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ABSTRACT 


The  potential  problem -is  solved  for  the  circulation  distribution  of  an 
optimum  propeller  with  a  finite  number  of  blades  and  a  hub  of  constant  diameter. 

The  effect  of  the  hub  has  been  calculated  for  specific  cases,  showing  that  it 
becomes  important  for  propellers  with  large  hub  diameters  and  small  number  of 
blades  and  increases  with  increasing  pitch. 

INTRODUCTION 

The  determination  of  the  bound  circulation  for  an  optimum  propeller  was  initially  per¬ 
formed  by  Goldstein1  who  was  able  to  relate  the  circulation  distribution  of  a  propeller  with  a 
finite  number  of  blades  to  one  having  an  infinite  number.  Goldstein’s  solution  was,  however, 
performed  for  the  case  of  a  propeller  having  a  zero  hub  diameter  and  when  the  vortex  sheet 
extends  to  the  propeller  axis.  Such  a  solution  is  considered  sufficiently  accurate  for  propellers 
with  a  hub  diameter  which  is  relatively  small  compared  with  the  propeller  diameter;  in  this 
case,  the  presence  of  the  hub  is  assumed  to  have  little  effect  on  the  circulation  distribution 
along  the  radius. 

However,  the  increasing  use  of  propellers  with  relatively  large  hubs  has  emphasized 
the  need  of  a  solution  for  such  cases.  The  problem  was  originally  investigated  by  Lerbs2 
who,  through  the  use  of  corrective  “ induction  factors,”  was  able  to  consider  the  presence 
of  an  infinitely  long  hub.  A  solution  to  the  potential  problem  has  also  been  given 
by  McCormick3  for  the  effect  of  a  finite  hub  but  with  simplified  boundary  conditions  given  at 
the  hub  radius.  The  effect  of  such  assumptions  is  to  produce  a  discontinuous  change  of  cir¬ 
culation  at  the  hub  which  is  contrary  to  available  experimental  evidence. 

The  problem  is  solved  for  a  propeller  having  a  minimum  energy  loss  and  for  which  the 
flow  far  behind  the  propeller  can  be  considered  to  be  the  same  as  that  formed  by  rigid  trailing 
vortices  moving  backwards  at  constant  angular  velocity.  The  hub  is  assumed  to  be  of  constant 
diameter  and  extending  from  the  propeller  plane  to  infinity. 

STATEMENT  OF  THE  PROBLEM 

The  problem  is  analogous  to  that  of  a  propeller  with  zero  hub1  but  with  modified  bound¬ 
ary  conditions.  In  this  case,  however,  solutions  which  have  singularities  at  the  shaft  axis 
can  also  be  considered. 


1R«ler«ncM  m  Hated  on  pago  16. 


Laplace’s  equation  for  this  type  of  flow  can  be  reduced  to 


*  +  (*  +  m2) 
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where  <f>  is  the  velocity  potential, 


u  is  equal  to  —  , 
v 

0  is  the  angular  velocity, 

v  is  the  velocity  of  advance, 

r  is  the  distance  from  axis  of  rotation, 

03 

£  is  equal  to  6  -  — ,  and 

v 

r,  0,  z  are  cylindrical  coordinates. 

The  boundary  condition  which  has  to  be  satisfied  on  the  vortex  sheets  is  that  the  velocity  be 
normal  to  the  sheet  and  can  be-  written  as 


d<f> 


on 


C  -  0  or  2lL 
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for  rx<r  <R  ot  <  fi  <  n0 


[2] 


where  p  is  the  number  of  blades, 
R  is  the  propeller  radius, 
rl  is  the  hub  radius, 

is  equal  to - ,  and 
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lil  is  equal  to- 
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Additional  conditions  have  to  be  imposed  at  r  -  R  where  continuity  should  exist  for  both  <f> 
and  grad  <f>.  Furthermore,  grad  <f>  should  also  vanish  at  an  infinitely  large  radius. 

The  presence  of  the  hub  imposes  certain  boundary  conditions  on  both  the  circulation 
at  the  vortex  sheet  and  on  the  induced  flow.  When  approaching  r1  from  inside  the  hub,  the 
circulation  has  to  be  zero.  When,  however,  rx  is  approached  from  outside  the  hub,  the  possi* 
bility  exists  that  the  circulation  may  change  either  discontinuously  or  continuously  to  zero. 

A  discontinuity  in  the  circulation  distribution  at  the  hub  will,  however,  mean  that  a  discon¬ 
tinuity  will  exist  in  the  pressure  field  between  two  adjacent  blades.  Physical  considerations 
preclude  the  possibility  of  the  existence  of  such  pressure  changes,  and  the  flow  will  tend  to 


2 
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equalize  from  the  pressure  side  of  one  blade  to  the  suction  side  of  the  other.  Furthermore, 
available  experimental  information4  indicates  that  the  radial  distribution  of  circulation  is 
continuous  to  zero  at  the  hub.  Hence,  the  boundary  condition  is  written  as 

</)  =  0  at  r  =  fj  or  ^  [3] 

The  presence  of  the  hub  also  requires  that  the  radial  component  of  the  velocity  be  zero.  By 
this  condition 


d<f> 

—  =0  at  fsf.  or 
<9/i 
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SOLUTION  OF  THE  POTENTIAL  PROBLEM 
THE  EXTERNAL  FIELD,  t>R 

The  conditions  of  the  problem  are  such  that  <£  has  to  be  an  odd  function  of  £,  single¬ 
valued  and  continuous.  Assuming  that  the  solution  can  be  expressed  in  a  sine  series,  we  find 
that  the  coefficient  of  sin  pn£  must  be  a  linear  function  of  / pn(Pnv)  anc*  K pnC?71#*)*  where  l  and 
K  are  the  modified  Bessel  functions.  But  / pn(Pnv)  cannot  occur  since  grad  <j>  must  vanish  when 
r,  or  (iy  is  infinite.  A  solution  independent  of  £  can  also  exist  which  satisfies  Equation  [1]  and 
which  is  in  the  form  of  log  Hence,  we  may  assume 


Kpn(PnH  o) 


sin  pn  £  +  c  log 


Po 


[5] 


where  cn  are  constants  to  be  determined  with  Kpn(pnp0)  and  log  pQ  being  inserted  in  the 
denominator  for  future  convenience. 


THE  INTERNAL  FIELD,  ^  <  r  <  * 

In  this  region  we  write 

a2 

<f>  ~  ~  - - 2  £  +  ^1  [6] 

1+p 

where  <£.  is  a  velocity  potential  to  be  determined.  For  such  a  case  the  boundary  condition 
given  by  Equation  [2],  becomes 

-**1  =  0  at  £  =  0  and 

K  P 

and  Laplace’s  equation 
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Since  £  is  periodic  in  the  interval  0  to  2  it /p,  it  can  be  expanded  in  a  half-range  cosine  series 


»  cos  (2  m  +  1)  - 
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and  for  the  same  interval,  <f>l  can  be  expanded  in  a  semi-infinite  cosine  series  in  £ 


<t>l=f0(p)+  £  fM  cos(2m  +  l)^ 

m  =  0 

Differentiating  term  by  term  and  substituting  in  Equation  [7] 
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and 


H)' ^  T^r)  “» 


Equation  [10]  can  be  integrated  to  give 
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where  D  log  n0  is  an  arbitrary  constant.  We  now  consider  a  new  function  such  that 


Ue>---£  -J— . 

*t  (2m +  1)2  |_l  +  <*  J 


[13] 


and  Equation  [11]  becomes 
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K) 


-(2m  +  l)2  (1+M2)-^-  gm(n)  +  /x 


2ll  ..2 
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[14] 


The  solution  of  this  homogeneous  differential  equation  can  be  recognized  as  a  Bessel  function, 
while  the  particular  solution  can  be  written*  as  a  Lommel  function.  The  general  solution  for 
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16p2(l-/x2) 


1  +  p2  p2(2m  +  l)2  (1+ n2)4 


[16] 


256  u2  (1  -14  fi2  +21  /x4  -  4  jt®) 

+  - - - + 


p4  (2m  +  l)4  (l  +  p2) 
bm'  =  arbitrary  constants 

Substituting  the  values  for  /0(fx),  fm(p)  into  Equation  [9]  and  hence  into  Equation  [6] 

we  get 


*-\z 


m  =  0 


4  T1)V(t^)  ^  Kv(vp) 


cos  v£  +  D  log —  [17] 


where  om  and  bm  are  arbitrary  constants  to  be  determined  from  the  boundary  conditions. 


DETERMINATION  OF  THE  ARBITRARY  CONSTANTS 

The  arbitrary  constants  om,  bm,  c  n ,  c,  and  D  are  determined  from  the  following  boundary 
conditions: 

(a)  is  continuous  at  r  =  R 
36 

(b)  -  is  continuous  at  r  =  R 

dp 

(c)  <f>  =  0  at  r  =  fj 
3<f> 


(<*) 


dfi 


0  at  r  =  t. 


In  order  to  evaluate  the  constants,  it  is  first  necessary  to  express  <£  for  the  external  field  in 
terms  of  sin  pn  £  in  the  interval  0  <  £  <  2n/p.  In  this  range  we  write 


cos  v 


<4  2 


j  4n2  -  (2m  +  l)2 


sin  pn£ 


[18] 
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Substituting  Equation  [18]  in  Equation  [17]  and  equating  to  Equation  [5]  for  p  =  p0  and  all  £, 
we  satisfy  boundary  condition  (a)  when 


8  \  _ n _  I" 4^  £i>Mq) 

^To  4n2-(2m  +  l)2  (2m  +  l)2  +  ^  +  mJ 

Similarly  for  boundary  condition  (b), 
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and 

Z?  =  c 

where  primes  denote  differentiation  with  respect  to  the  argument. 
For  zero  circulation  at  the  hub 
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and  from  (d) 
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Equations  [19],  [20],  [21],  and  [22]  can  be  solved  for  the  coefficients  am,  bm,  c„,  and  D.  From 
Equation  [19]  and  [20],  we  obtain 
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The  determination  of  the  constants  am  and  b  is  now  considered  from  Equations  [21]  and 
[23].  A  method  of  successive  approximations  can  be  used  to  determine  the  coefficients  from 
this  infinite  set  of  simultaneous  equations.  However,  the  singularity  which  exists  at  the  edges 
£  =  0  and  2  n/p  will  result  in  a  very  slow  convergence.  Restricting  the  solution  for  values  of 
fiQ  which  are  not  too  small,  we  can  write  from  Equation  [16] 


Ti,  1/(^0)* 


This  approximation  improves  with  increasing  number  of  blades.  Also 
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Furthermore,  using  Nicholson’s6  asymptotic  expansions 


Kpn  (Pn M0)  (1+Mo2)* _ Mo _ Eq(*-/*o2) 

Kpn(PnH)  H  2?n(1+#i o2)  8P2"2  (l  +  fto2^2 

or  for  /z0  not  too  small 

Kp'n(PnPo)  ~  (l-t-^02)^ 

Kpn(PnP  0)'”  PO 

Similarly 

V(^o)  (1  +  Po>*  Po  M4~Po2) 

-  -n  —  "  +  5/2  “ . 

^(‘Vo)  Po  ^(l+Mo2)  8*/(1+f‘o2) 


_  (1+Eq2)* 
Po 

Equation  [23]  can  then  be  written 
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For  usual  values  of  fiQ,  it  is  possible  to  neglect  the  second  term  on  the  right-hand  side  of 
[23]  as  compared  with  the  first.  Then, 


2T — — 
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and  from  Equation  [21] 
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and  finally  we  obtain 
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for  n  =  l,  2,  3,  .  .  . 

The  existence  of  a  solution  for  this  infinite  system  cannot  be  shown,  as  the  necessary  theorems 
have  conditions  which  are  too  restrictive  in  this  case.  However,  the  system  is  convergent,  and 
a  solution  can  be  found  by  the  standard  method  of  successive  approximations. 

It  will  be  noted  that  as  the  hub  diameter  tends  to  zero,  i.e.,  nl  tends  to  zero,  Equation 
[27]  degenerates  into  a  system  of  equations  given  by  Goldstein. 


DETERMINATION  OF  THE  BOUND  CIRCULATION 


The  circulation  distribution  is  proportional  to  the  discontinuity  in  velocity  potential  at 
the  screw  surface  {  =  0  or  2n/p.  If  [0]  denotes  the  discontinuity  in  the  potential,  then  from 
Equation  [17] 

[*]=£ 

p 

It  is  then  possible  to  write 
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If  T",  is  the  circulation  of  a  blade  section  for  a  propeller  having  an  infinite  number  of  blades, 
then 


vj*£oo 
2  ntov 
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Defining 
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is  the  factor  for  a  propeller  having  a  finite  hub.. 
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or  from  Equation  [26] 
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EFFECT  OF  FINITE  HUB 


The  use  of  Equation  [27]  for  the  determination  of  the  coefficients  am  involves  consider¬ 
able  labor,  and  in  order  to  indicate  the  effect  of  the  hub,  an  alternative  set  of  equations  has 
been  derived  in  the  Appendix.  These  equations  involve  simplifications  of  the  basic  development 
but  have  good  accuracy  for  the  usual  values  of  nQ  and  for  hub  diameters  not  larger  than  half 
the  propeller  diameter.  The  errors  which  are  introduced  decrease  considerably  with  increasing 
number  of  blades. 

Through  the  use  of  Fquations  [31],  [39],  [40],  and  [43]  of  the  Appendix,  calculations  were 
performed  showing  the  effect  of  hub  diameters  of  0.2,  0.3  and  0.4  for  a  four-bladed  propeller. 

Both  the  factor  *  and  the  circulation  distribution  are  shown  on  Figures  1  and  2  for  this  propeller 
operating  at  1/X.;  =  ft 0  =  4.0.  Note  that  the  circulation  function  falls  off  to  zero  at  the  hub  cor¬ 
responding  to  the  boundary  condition  (c)  for  the  circulation.  Near  the  tip,  the  distributions 
for  finite  and  zero  hub  coincide;  this  is  to  be  expected  since  the  influence  of  the  hub  must 
decrease  towards  the  tip.  Figure  3  shows  the  ratio  of  the  factor  k  for  a  finite  hub  and  a  zero 
hub  for  various  hub  diameters. 

The  effect  of  number  of  blades  on  k  is  shown  on  Figure  4  for  a  two-,  a  four-and  a  six- 
bladed  propeller  operating  at  1/A  .  =  /i0  =  4.0  for  r/R  =  0.2.  Note  that  the  effect  diminishes 
with  increasing  number  of  blades. 

It  should  be  noted  that  although  the  presence  of  the  hub  has  an  appreciable  effect  on 
the  circulation  near  the  hub,  the  velocities  are  the  same  over  the  entire  radius.  Both  the 
condition  of  normality  and  Betz’s  condition  hold  for  a  finite  hub,  and  they  are  sufficient  to 
uniquely  determine  the  induced  velocity  components. 
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p«r 


Figure  2  •  Effect 


Figure  3  •  Effect  of  Hub  Diameter  for  *„/ Kh  for  nQ  =  4,  p  =  4 


Figure  4  *  Effect  of  Number  of  Blades  for  nQ  ■  4,  r^/ft  -  0.2 
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CONCLUSION 


The  solution  of  the  potential  problem  for  a  propeller  with  finite  hub  indicates  that  the 
presence  of  the  hub  results  in  a  decreased  circulation  at  the  inner  radii.  This  effect  becomes 
important  for  propellers  with  large  hub  diameters  or  small  number  of  blades  and  increases  with 
increasing  pitch. 


APPENDIX 

APPROXIMATE  SOLUTION  FOR  THE  COEFFICIENTS 


The  solution  of  Equation  [27]  for  the  coefficient  am  involves  considerable  labor,  and 
an  attempt  has  been  made  to  approximate  the  solution  of  this  system  for  usual  values  of 
For  this  purpose,  let 


Mo 


am  =  - 

tn 


l^o2  m 


A_  +  8  A, 


[30] 


where  Am  is  the  coefficient  defined  by  Goldstein1  and  is  given  by 

(2m)  ! 


22m  (m! )2  (2m  +  1) 


[31] 


and  SA  is  a  correction  coefficient  to  be  determined.  Furthermore,  Nicholson’s  asymptotic 
expansions6  are  used  for  the  /  and  K  functions,  where 


and 


where 


[32] 


[33] 


[34] 


Substitution  of  Equations  [30],  [32],  and  [33]  into  the  system  of  Equation  [27],  gives 


12 


Zu. 

m  =  0 


r 


i 

-(2m  +  l) 


2n  +  (2m  +  l)_  ir  \l  +  #*02/ 


OO 

27 

m  =  0 


-v(,y0-y,)  m 
e  0  1  r, 


[2n+(2m  +  l)]  (2m  +  l)T‘ 


27 


4 


-2v(y0-y1) 

,e 

2n  +  (2m  +  1) 


[35] 


For  usual  values  of  fiQ  and  for  values  of  r^/R  which  are  smaller  than  0.5,  Am  e  0  1  is 

small  compared  with  Tl  v(v n^)/ (2m  +  l)2.  Furthermore,  the  second  term  of  the  left-hand  side 
of  Equation  [35]  is  also  small  compared  with  the  first  and  can  safely  be  neglected.  Equation 
[35]  can  then  be  written 


z 


M. 

2n-(2m  +  l) 


Wl  +  ^\*  y  t  V(y°  yi) 

[2n  +  (2m  +  l)]  (2m +  1)2 


[36] 


Using  Goldstein’s  definition  of  the  T  function  as 


f*  l 


-  F_ 


1  +Mi2  "2’  (2 »+  4) 


[37] 


it  is  found  that  the  terms  corresponding  to  values  of  m  =  1,  2,  3  ...  of  the  right  hand-side  of 
equation  can  be  omitted  with  good  accuracy.  The  error  introduced  by  such  an  approximation 
is  less  than  3-percent  ir.  the  determination  of  the  value  of  Am,  for  a  two-bladed  propeller, 
rx/R  =  0.4  and  fiQ  =  3.  The  error  decreases  for  increasing  number  of  blades  and  increasing 
Equation  [36]  becomes 


27 


for  n  =  1,2,3  .  .  . 


n."a<y°’yi> 

(2n  +  l) 

[38] 


The  F  functions  for  two  and  four  blades  are  plotted  on  Figures  5  and  6.  Let 


■-CjHiV t'H 


8A„  = 


s  2  B, 


[39] 


then  B  m  has  to  satisfy  the  equation 


27 


m  »  0 


2n-(2m  +  l)  2n  +  l 


for  n  =  1,2,3.  .  . 


[40] 
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Solution  of  the  system  [40]  gives 


B0  =  0.4762 
Bx  =  0.1128 
B2  =  0.0531 

53  =  0.0309 

54  =  0.0200 


Bs  =  0.0136 
B6  =  0.0095 
B7  =  0.0066 
Bs  =  0.0044 
Bg  =  0.0025 


When  am  is  represented  by  Equation  [30],  the  equation  for  the  factor  ka,  becomes 


ri,»>P)  Vp  A 

(2ot  +  1)2  1  +  Mq  ^v^vVq)  _ 


m  =  0  L 


liiE 


(^i) 


(2m  +  l)2 


M°  4  tvft) 

Kv(vnx) 


1  +  ■• 2 


[41] 


2  A+M*\  5^  S/4 

n\  V2  •  7?o  LM^o)  Kv{vh)  Iv{vh) 


It  is  noticed  that  the  first  term  of  Equation  [41]  is  equal  to  the  ratio  of  the  circulations  for  a 
propeller  having  zero  hub.  Denoting  this  by  kq  and  writing 


“A*** o  " 


then 


5k 


2,/i v\  y  r 4  4  mo  /, 

M  Ln  (2w*  + 1)2  mi+Mo2 

"L/i>**o)  m*vi>  /v(v^0)j  r 


y(yMi)  I 

(*V  J 


[42] 


[43] 


The  factor  5k  can,  therefore,  be  considered  as  a  correction  factor  to  the  case  of  a  propeller 
with  zero  hub  and  can  be  evaluated  using  Equations  [31],  [39]  and  [40]  for  the  coefficients 
Am  and  54  „. 

in  m 
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;tion  F  for  p  ■  4 
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